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ABSTRACT 

The effect of the stellar flux on exoplanetary systems is becoming an increasingly important property 
as more planets are discovered in the Habitable Zone (HZ). The Kepler mission has recently uncovered 
circumbinary planets with relatively complex HZs due to the combined flux from the binary host stars. 
Here we derive HZ boundaries for circumbinary systems and show their dependence on the stellar 
masses, separation, and time while accounting for binary orbital motion and the orbit of the planet. 
We include stability regimes for planetary orbits in binary systems with respect to the HZ. These 
methods are applied to several of the known circumbinary planetary systems such as Kepler-16, 34, 
35, and 47. We also quantitatively show the circumstances under which single-star approximations 
break down for HZ calculations. 
Subject headings: astrobiology - planetary systems 
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1. INTRODUCTION 

The Habitable Zone (HZ) is broadly denned as the 
region around a star where water can exist in a liquid 
state on the surface of a planet with sufficient atmo- 
spheric pressure. The boundaries of the HZ are usually 
calculated based upon the properties of the host star, 
with assumptions regarding the response of the plane- 
tary atmosphere to stellar flux. Attempts to quantify 
these boundaries have evolv ed cons i derably sinc e i nitial 
discus sions on the topic by iHuangj ()1959l . I1960D . IHartl 
(1979) provided estimates for the Solar System and later 
models of runaway greenho use for Venus w ere incorpo- 
rated int o the calculations bv | P ollack (197j]) and IKasting! 
(11988ft . IKasting et all 11993J) calculated detailed one- 
dimensional (altitude) climate models and considered 
conditions whereby the equilibrium would sway to a run- 
away greenhouse effect or to a runaway snowball effect, 
thereby providing robust estimates of the HZ boundaries 
for main sequence stars. 

Interestingly, much of the work on the HZ preceded 
the discovery of exoplanets. The exoplanet parameter- 
space, specifically orbital period and planetary mass, 
is slowly being expanded by different techniques. For 
example, the radial velocity (RV) technique has ex- 
plored well beyond the HZ for main sequence stars, 
allowing analysis of H Z properties of these systems 
(jKane fc Gelinoll2012aD . Exoplanet systems tend to be 
dominated by Jovian planets, many of which are in 
highly eccentric orbits which causes t he planet to pass 
through the HZ (jKane fc Gelinoll2012bT ). There have also 
been several notable d iscoveries of supe r-Earths in the 
HZ, such as Gl 581d (llldrv et al.l 12007^ an d Gl 667Cc 
([Anglada-Escude et al.ll2012t iDelfossd 120121 ). The tran- 
sit method is also contributing to the census of HZ plan- 
ets. The Kepler mission is de signed to explore the HZ 
(|Kaltenegger fc Sasselovl 1201 If) and has already found 
several which meet this criteria, such as Kepler-20b 
(|Borucki et al.|[2012T) . 

One of the more interesting discoveries to result from 
the Kepler mission is that of circumbinary planets, or 
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planets that orbit the center-of-mass of a central binary 
system. The first of these announced was Kepler- 16b 
([Doyle et al.H201lD. soon to be followed by Kepler-34b & 
Kepler-35b ([Welsh et al.ll2012ft. and t he multi-planet sys- 
tem Kepler-47b.c (jOrosz et aL| [2~012'). There have been 
several discussions regarding the HZ for these systems, 
though these to-date approximate the calculations as a 
single -star for the purposes of determining th e HZ bound- 
aries (jQuarles et al.ll20"rl lOrosz et al.ll20rl . 

Here we perform a more thorough analysis of the HZ 
boundaries for circumbinary planetary systems. We con- 
sider the contributions of both stars and the resulting 
position dependence of the HZ boundaries. We thus de- 
termine the mass/separation ratio dependence of these 
boundaries and include the orbital motion of the binary 
to produce a time dependent map of the HZ. We ap- 
ply these calculations to the known Kepler circumbinary 
planetary systems of Kepler-16, 34, 35, and 47. 

2. HABITABLE ZONES OF SINGLE STAR SYSTEMS 

The HZ for single star systems defined by stellar irra- 
diation has been discussed by a number of sources, most 
notably by IKasting et al.l ([1993). The boundaries of the 
HZ for main sequence stars have been further quantified 
as a function of effective temperature ( Underwood et al.l 
[2001 [Se¥is~eFaT] [20071: bones fc SleeplboiOft . Note that 
there are additional in fluences on the hab i tability such 
as th e effects of tides (IBarnes et al.l [20091: IHeller et al.l 
120 111) and orbital stability ([Kopparapu fc Barnesll2010ir 
Here we consider only the irradiation effects on the HZ 
and briefly describe the essential components which we 
expand upon in future sections. 

One of the fundamental stellar properties which deter- 
mines the extent of the HZ is the luminosity of the host 
star, which is approximated as 



off 



(1) 

where a is the Stefan-Boltzman n constant. The model 
adopted by IKasting et al.l ()1993D utilizes the principle of 
a planetary energy balance: net incoming solar radia- 
tion must equal the net outgoing infra-red (IR) radiation. 
This method determines the stellar flux, S, for which 



2 



Stephen R. Kane & Natalie R. Hinkel 



the outgoing IR radiation is sustainable for various at- 
mospheric assumptions. The underlying implication here 
is that the stellar flux received at the HZ boundaries is 
dependent on the amount of IR radiation incident upon 
the upper atmosphere. Thus there is a temperature de- 
pendence on the host star. Using the boundary condi- 
tions of runaway greenhouse and maximum greenhouse 
effects at the inner and outer edges of the HZ respec- 
tively (jUnderwood et al]|2003f ). the stellar flux at these 
boundaries are given by 



5 inner = 4.190 x l(T 8 T c 2 ff - 2.139 x l(T 4 T f 



'off + 1.268 (2) 

Pouter = 6.190xl0- 9 T c 2 ff -1.319xl0" 5 T cff +0.2341 . (3) 

The inner and outer edges of the HZ are then derived 
from the following 
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where the radii are in units of AU and the stellar lu- 
minosities are in solar units. The dependence of these 
boundaries upon the wavelength of the incident radiation 
is of critical importance when considering two separate 
radiation sources. 

3. HABITABLE ZONES OF MULTI-STELLAR SYSTEMS 

In this section we generalize the above principles to 
multi-stellar systems. The HZ of planets in S-type orbits 
(planets orbiting one compon ent of a ste l lar bi nary) has 
been considered in detail by lEggl et all (|2012D and ap- 
plied to a Centauri B by iForganl (12012ft. now know n to 
host a terrestrial exoplanet (Dumusq ue et al.ll2012D . We 
address the specific case of HZ boundaries for planets in 
P-type orbits, where the planet orbits the center of mass 
of both stars. 

3.1. Combined Black-body Spectrum and Stellar Flux 

The radiation received at any point in a binary system 
is a combination of two blackbody spectral radiances. 
First we consider Planck's law as a function of wave- 
length, A: 

J (A,D = ^ ehe/ J r _ 1 (6) 

where T is the effective temperature. Note that this 
equation is per steradian. The stellar flux received at 
a particular location may then be expressed by the fol- 
lowing equations for star 1 and star 2: 
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where R+ i and i?*,2 are stellar radii and r\ and r 2 
are the distances from a given location to star 1 and 
star 2 respectively. The total flux received at this lo- 
cation is then S — Si + S 2 - Wien's Displacement Law 
(A max 7cff = 2.9 x 10 -3 ) may then be used to estimate the 
equivalent effective temperature of a single energy source 
that would produce the same energy flux. For the pur- 
poses of integrating the total flux over all wavelengths, 
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Fig. 1. — The blackbody spectra of the individual stars in the 
binary system and the combined flux for both stars. The peak of 
the combined spectrum peaks at a wavelength of 528 nm (vertical 
line) which, using Wien's Displacement Law, is equivalent to a 
temperature of 5481 K. 

Equations [7] and [5] reduce to the Stefan-Boltzmann law. 
However, it is important to consider the wavelength de- 
pendence of the incident flux, as we will see in later sec- 
tions. 

As an example, consider the case of a binary system 
consisting of a G2V star and a K5V star, separated by 
0.1 AU. Shown in Figure Q] are the individual and com- 
bined blackbody spectra for this system. The total flux 
peaks at a wavelength of 528 nm which is equivalent to a 
temperature of 5481 K. This particular evaluation of the 
blackbody spectra does not take into account the dis- 
tances to the stars. When the radii and distances are 
included, this introduces a spatial dependence. 

Figure [2] shows a flux map of the same system with 
two different zoom levels to emphasize the contribution 
from both stars. For the purposes of contrast, the calcu- 
lated flux has been converted to a logarithmic scale and 
normalized to lie between and 1. On this flux scale, 
the contours are equal to values of 0.1 to 0.5 in steps of 
0.1 moving from the outer edge inwards. There is a clear 
asymmetry associated with the presence of the second 
star which is particularly noticeable at distances closer 
than 0.3 AU. 

3.2. Habitable Zone Boundaries 

The traditional HZ due to stellar irradiance for a sin- 
gle star is calculated as a function of both the stellar flux 
received and the peak wavelength of the energy distribu- 
tion. For a single star, this has a radial dependence in 
a spherical geometry which allows an analytical solution 
for the inner and outer HZ boundaries. For a multi- 
stellar system, the locations of these boundaries is far 
more complex, since both the stellar flux and peak wave- 
length have non-radial dependencies due to the combina- 
tional effect produced from the two stars at a particular 
location. 

We approach this problem by constructing a two- 
dimensional grid in the orbital plane of the binary star. 
For each grid location we perform the following steps: (1) 
calculate the combined stellar flux, (2) calculate the com- 
bined blackbody function and determine the peak, (3) 
convert the wavelength for the peak into effective tem- 
perature using Wien's Displacement Law, and (4) calcu- 



3 




Fig. 2. — The stellar flux map, depicted with both gray scale and contours, of a binary system consisting of a G2V and a K5V star with 
a separation of 0.1 AU. The right panel is a zoom of the left panel. In each panel, the calculated flux has been converted to a logarithmic 
scale and normalized to lie between and 1. On this scale, the contours are equal to values of 0.1 to 0.5 in steps of 0.1 moving from the 
outer edge inwards. 
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Fig. 3. — The stellar flux map shown as a gray scale map (as for Figure[2} and the HZ bound aries for two different binary configurations. 
The left panel depicts these quantities for the G2V-K5V binary discussed in detail in Section l3.ll where only the inner edge of the HZ is 
shown on this scale. The right panel depicts these for an equal-mass binary consisting of two MOV stars separated by 0.5 AU. In this case 
both the inner and outer edge of the HZ are shown on the plot with a very clear asymmetry to those boundaries. 



late the expected stellar flux at the HZ boundaries for 
this temperature using Equations [2] and [3] If the calcu- 
lated stellar flux at that location matches that expected 
for a HZ boundary then that grid position is marked as 
such. 

Figure [3] depicts the HZ boundaries for two example 
binaries: the G2V-K5V binary described in ij3.ll with a 
separation of 0.1 AU (left panel), and an equal-mass bi- 
nary consisting of two MOV stars separated by 0.5 AU 
(right panel). The inner HZ boundary for the G2V-K5V 



binary appears symmetric on this scale, but the distance 
of the boundary from the primary (G2V) star varies be- 
tween 0.90 and 0.94 AU. For the secondary star, this dis- 
tance ranges from 0.84 to 1.01 AU. The asymmetric HZ 
boundaries for the M0V-M0V binary shown in the right 
panel represents a far more extreme case where both the 
inner and outer boundaries are plotted. 

The asymmetry in the HZ boundaries may be quanti- 
fied as a function of the binary separation by considering 
the axis ratio of both the inner and outer boundaries of 
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Fig. 4. — The dependence of the HZ asymmetry on the binary separation for the previously discussed cases of a G2V-K5V and MOV— MOV 
binary. At each separation, the axis ratio of both the inner and outer HZ boundaries are computed. 



the HZ. We performed a simulation of this dependence 
by calculating the HZ boundaries in the cases of the two 
examples described above. We utilized a range of separa- 
tions between 0.01 and 1.0 AU and located the minimum 
and maximum values of both the inner and outer bound- 
aries. The results of this simulation are shown in Figure 
[H The case of the G2V-K5V shows moderate variation 
in the inner HZ asymmetry but negligible variation in 
the outer HZ. For the equal-mass M0V-M0V binary, the 
asymmetry rapidly increases with increasing separation 
for both the inner and outer boundaries. 

3.3. The Keplerian Orbit of the Binary 

The orbital elements of binary stars are known to 
encompass a broad range of orbital perio d s and mass 
ratios (see fo r exam ple lHalbwachs et al.l (|2005() and 
iMathieu et~aT1 (|2004D V This diversity results m an 
equally varying range of HZ structures due to the mass 
ratios and binary separations. From a fixed location 
relative to the binary center-of-mass, the apparent HZ 
boundaries will have a time dependence owing to the or- 
bital motion of the binary. 

For a Keplerian (non-circular) binary orbit, assuming 
a fixed reference point, the changing HZ boundaries will 
have an additional disparity resulting from the varying 
binary separation. This will result in an oscillating HZ 
structure as the binary rotates around their center-of- 
mass. As shown in the previous sections, the oscillation 
will have a much larger effect for low-mass stars than 
for solar-type stars. The large variation in flux within 
typical HZ regions may ultimately render some binary 
systems uninhabitable. 

4. PLANETS WITHIN THE HABITABLE ZONE 

Now we consider a planetary orbit which is 
passing through asymmetric HZ regions and the 
flux/temperature variations that result. 

4.1. Stable Configurations 

Stable planetary orbits in both S-type and P-type sys- 
tem s have been discus sed by a variety of authors, includ- 
ing iHarringtonl (Il977l). but more r ecent ly by such au- 
thors as lEggleton fc Kiseleval (|1995t ) and Musiela k et al.l 




Fig. 5. — A contour plot of the dependence of the critical semi- 
major axis of a planetary or bit, q c , on the mass-ratio , fi, and sep- 
aration, aj, of a binary star (Holman & Wiegcrt 1999). 

(2005). We adopt the formalism of IHolman fc Wiegerti 
(1999) which provides analytical solutions for planetary 
orbital stability in binary star systems. For P-type or- 
bits, this is a function of the binary separation, a^, the 
mass ratio, /i, and the eccentricity of the binary orbit. 
The mass ratio is defined as /x = 1712/ (mi + 1TI2), such 
that an equal mass binary has a mass ratio of [i = 0.5. 
For a binary in a circular orbit, or e = 0, Equation 3 of 
IHolman fc Wiegerti (jl999) may be expressed as 

a c = (1.60 + 4.12 j u-5.09/i 2 )a b (9) 

where a c is the minimum allowed semi-major axis of the 
planet. Figure [5] shows the dependence on a& and fj, 
of the critical semi-major axis of a planet orbiting the 
center-of-mass of the binary. For binary separations of 
ab < 0.2 AU, the critical semi-major axis is only weakly 
depend ent upon the mass ratio of the binary. As pointed 
out by IHolman fc Wiegerti (|1999f) , there often exists in- 
stability islands beyond a c which occur at mean-motion 
resonances. Thus the critical semi-major axis calcula- 
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Fig. 6. — The stellar flux map shown as a gray scale map and HZ boundaries (solid lines) for a G2V-K5V (left panel) and MOV-MOV 
(right panel) binary, as for Figure [3] Also shown in each case is the boundary for the critical semi-major axis (dotted line), beyond which 
planetary orbits may retain long-term stability. This shows that the stable planetary orbits may exist in the HZ for the G2V— K5V binary, 
but stable orbits are excluded for the HZ of the MOV-MOV binary. 



tions here represent a lower limit on the total stability 
regions within a given system. 

Applying this stability criteria to the two examples dis- 
cussed thus far will yield whether or not a stable plan- 
etary orbit can be maintained within the HZ of those 
systems. Figure [5] is similar to Figure [3] in that it shows 
the stellar flux map and HZ boundaries for the G2V-K5V 
and MOV-MOV binaries, but slightly zoomed out so the 
outer part of the HZ for the G2V-K5V binary is visible. 
Also shown are dotted lines which indicate the critical 
semi-major axis for planets within the system, centered 
upon the binary center-of-mass in each case. This clearly 
shows that stable planetary orbits within the HZ for the 
G2V-K5V binary are possible but the the presence of 
planets are excluded from the HZ of the MOV-MOV bi- 
nary. 

4.2. The Keplerian Orbit of the Planet 

A planet in orbit around the center-of-mass of the bi- 
nary will experience variable conditions due to the rota- 
tion of the binary, as mentioned in i )3.3l An additional 
component of the total flux variation is due to the Kep- 
lerian motion of the planetary orbit. The radial depen- 
dence of the flux means that the influence of the binary 
motion received at the planet will be more pronounced 
when the planet is passing through periastron passage. 

As a demonstration, we include the presence of a planet 
in the G2V-K5V binary system example used through- 
out this paper (recall that our particular MOV-MOV bi- 
nary example excludes stable orbits within the HZ for 
the binary separation used). We show the stellar flux 
map (see §3.1[) in Figure along with the HZ bound- 
aries (see §3.2[) . the critical semi-major axis boundary 
(see N4.ip . and the orbit of the planet. The planet has 
a semi-major axis of 0.9 AU and an eccentricity of 0.6 
such that it spends most of the orbit within the HZ of 
the system. At each location of the orbit, we recalculate 



the flux received by the planet and estimate the equilib- 
rium temperature of the planet assuming that the inci- 
dent flux is red i stribu ted around the entire planet (see 
iKane fc Gelinol ()2011l )). We also account for the stel- 
lar binary motion at each location during the planetary 
orbit. The bottom panel of Figure [Jj shows the change 
in this temperature during one complete orbital phase 
as the planet moves from apastron (phase = 0.0) to pe- 
riastron (phase = 0.5) and back again. The long-term 
variation is caused by the Keplerian orbit of the planet 
and the short-term variations are caused by the 9 day 
orbital period of the binary. The temperature variations 
can be dominated by either the orbits of the binary or the 
planet, depending on the relative orbital configurations. 

4.3. Binary/Planet Dynamical Interactions 

Thus far we have assumed a decoupling of the Keplar- 
ian motion of the binary and the Keplerian motion of 
the planet. However, the secular dynamical interactions 
in the system will result in a break-down of the periodic 
behaviour one may otherwise obtain when you assume 
the planet and binary are orbitally decoupled. Although 
the semi- major axis of the planet's orbit will remain rel- 
atively unchanged, the shape (eccentricity) and orienta- 
tion (argument of periastron) will have time-dependent 
components. This has been described in detail by var- 
ious au thors, including analytical solutions fo r stellar 
triples ([Soderhielml 11984 iGeorgakarakosI 120091) . single- 
star multi-planet systems dLee fc PealeH2003|). effects o f 
general relativity (lEggleton fc Kisseleva-Eggletoni r2006) , 
and the dy namics and stability of p lanets in circumbi- 
nary orbits (poplin fc Blundelfeoill ). The amplitude of 
these variations depends on numerous factors, such as 
the r elative inclinations of the orbits (Fara gcTfe Laskarl 
[2011 . Here we assume co-planar orbits, as is the case 
for the Kepler systems discussed in the following section. 

The timescale of the secular dynamical interactions 
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Fig. 7. — Top panel: The G2V-K5V binary system showing the 
HZ boundaries (solid lines), critical semi-major boundary (dotted 
line), and a planetary orbit (dashed line). The planetary orbit has 
a semi-major axis of a = 0.9 AU and an eccentricity of e = 0.6. 
Bottom panel: The equilibrium temperature of the planet as a 
function of planetary orbital phase. Phase 0.0 is at apastron and 
phase 0.5 is at periastron. 

is largely dominated by the periastron precession of 
the planet which occurs on much shorter timescales 
than the precessio n of th e binary. The simulations by 
iDoolin fc BlundelH (|2011l ) of test particles in circumbi- 
nary orbits show that the period of the periastron pre- 
cession of an orbiting planet is related to the semi-major 
axis of the orbit, a, via the power-law oc a 3 ' 5 . Therefore, 
one can expect to be able to measure the effects of peri- 
astron precession for systems similar to those described 
here over the course of years, with a complete periastron 
precession period in the range of several decades. Al- 
though we are primarily concerned with the HZ of the 
binary systems, the stability and secular dynamics of the 
planetary orbits in those systems is an effect that one 
should consider when predicting the percentage of the 
orbital phase within the HZ over long periods or time. 

5. APPLICATION TO KNOWN SYSTEMS 

Here we apply the HZ boundary methods to several 
of the circumbinary exoplanetary systems found via the 
Kepler mission. 

5.1. Kepler-16 

Kepler- 16b was discovered by iDovle et al.l (|2011[) and 
was the first circumbinary planet announced as a result 
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Fig. 8. — Top panel: The Kepler-16 system showing the HZ 
boundaries (solid lines), critical semi- major boundary (dotted line), 
and planetary orbit (dashed line). Bottom panel: The equilibrium 
temperature of the planet as a function of planetary orbital phase. 
Phase 0.0 is at apastron and phase 0.5 is at periastron. 

of transit observations by the Kepler mission. The binary 
system consists of 0.7 and 0.2 Mq components in a 41 day 
eccentric orbit. The binary pair are orbited by a 0.3 Mj 
planet with an orbital period of 229 days. The system 
has been the subject o f numerous follow-u p studies, such 
as that carried out by iWinn et al.l (|2011l ) which showed 
that the plane of the stellar orbit, the planetary orbit, 
and the primarys rotati on are all closely alig ned. Radial 
velocity observations bv lBender et al.l (|2012| ) allowed the 
dynamical determination of masses for both of the binary 
components. 

Figure [5] shows the flux map, HZ boundaries, critical 
semi-major axis (stability) boundary, and planetary orbit 
for Kepler-16. The effective temper ature of the primar y 
star, T cS .i = 4450 K, is provided bv IDovle et "alt (|2011| ). 
For the secondary star, we adopt the te mperature of 
the sp ectral template star proxy used by I Bender et al.l 
(|2012l) . GJ 905, of T effj2 = 2800 K. Although there is a 
slight asymmetry in the inner HZ boundary, this is in- 
side the stability boundary. Overall, the flux is clearly 
dominated by the primary star which produces an outer 
HZ boundary centered on the primary. Our simulation 
shows that the offset between the center-of-mass of the 
system and the flux distribution results in the orbit of the 
planet moving in and out of the HZ during one complete 
orbit. In other words, the asymmetry in the temperature 
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Fig. 9. — Top panel: The Kepler-34 system showing the HZ 
boundaries (solid lines), critical semi-major boundary (dotted line), 
and planetary orbit (dashed line). Bottom panel: The equilibrium 
temperature of the planet as a function of planetary orbital phase. 
Phase 0.0 is at apastron and phase 0.5 is at periastron. 

variations experienced by the planet can occur even for 
circular orbits (see lower panel of Figure [5]) . 

5.2. Kepler-34 & Kepler-35 

IWelsh et al.l (|2012f) announced the discovery of two cir- 
cumbinary planetary systems from the Kepler mission: 
Kepler-34 and Kepler-35. These systems arc quite com- 
parable in terms of the mass ratio and separation of the 
binary stars. Figures l9l and ITOl show our calculations for 
the HZ conditions in both systems, respectively. The 
similarity of the stellar components in both cases yield 
analogous centroids for the HZ and stability boundaries, 
unlike the case of Kepler-16. The Kepler-34 binary con- 
sists of two near-solar components with T e g^ = 5913 K 
and T c ff 2 = 5867 K separated by 0.23 AU. Only the inner 
HZ boundary is shown since the combined flux of these 
stars pushes the outer boundary much farther out than 
for our own Solar System. The planet is in an eccentric 
orbit inside of the inner HZ boundary with a semi-major 
axis of 1.09 AU. The resulting planetary temperature 
variation has both clear short-term and long-term com- 
ponents. 

By comparison, the Kepler-35 system consists of stel- 
lar components with effective temperatures of 5606 K 
and 5202 K, respectively, separated by 0.18 AU. The 
HZ boundaries are thus smaller in size compared with 
Kepler-34. The planet's near-circular orbit is very close 
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Fig. 10. — Top panel: The Kepler-35 system showing the HZ 
boundaries (solid lines), critical semi-major boundary (dotted line), 
and planetary orbit (dashed line). Bottom panel: The equilibrium 
temperature of the planet as a function of planetary orbital phase. 
Phase 0.0 is at apastron and phase 0.5 is at periastron. 

to the stability boundary of the system, therefore the bi- 
nary. This causes the temperature profile of the planet 
to be very sensitive to the orbital motion of the binary, 
as can be seen in the plot of the temperature variations 
during one complete orbital period. 

5.3. Kepler-47 

The Kepler-47 system, discovered by lOrosz et al.l 
(2012), contains a binary star whose components are 
quite different, similar to Kepler-16. The two stars have 
temperatures of 5636 K and 3357 K, respectively, and 
are separated by only 0.08 AU. The resulting stellar flux 
dominated by the primary produces asymmetry in the 
flux only in the vicinity of the binary (within the sta- 
bility boundary). The main difference between Kepler- 
47 the other studied Kepler systems, however, is that 
this system has two known transiting planets. Figure [TT] 
shows our calculations for the Kepler-47 system, as well 
as the orbits and temperature profiles for the two planets. 
Both planets are in near-circular orbits with semi-major 
ax es as 0.30 AU and 0.99 AU, respectively. As noted 
bv lOrosz et al.l (j2012j ). the outer planet resides within 
the HZ of the system throughout the entire orbit. Since 
the planets are in non-eccentric orbits, the temperature 
variations are dominated by the binary orbit, with an 
additional variation caused by the orbital center-of-mass 
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offset from the primary stellar flux. The amplitude of 
the variations for the outer planet are relatively small, 
~ 5 K, which is negligible for habitability considerations 
in this region. 



6. CONCLUSIONS 

Calculating the amount of time an exoplanet spends 
within the HZ around a circumbinary system is rela- 
tively complicated. First, we must consider the amount 
of stellar flux from both stars received by the exoplanet 
at any point in its orbit fi )3.ip . Then we must determine 
the wavelength at which this combined flux peaks and 
how that translates into an effective temperature (£ I3.1|) . 
This will then allow us to determine the inner and outer 
boundaries of the HZ f b!3.2j) . Depending on the configura- 
tion of the binary masses and separation, the HZ bound- 
aries may then experience a significant time-dependent 
orbital motion with respect to the rotation of the stars 
( §3.3p . Finally, the stability of the binary system must 
be analyzed ( §4.1|) in order to fully understand the plan- 
etary orbit with respect to the oscillating HZ ( §4.2(1 . 

We have applied our analysis to a number of known 
exoplanet- hosting binary systems ($5]), namely Kepler- 
16, 34, 35, and 47. We have demonstrated that the cir- 
cumstances under which a single-star approximation is 
viable is solely reliant on the stellar masses and sepa- 
ration of the binary. For example, the components in 
Kepler-16 and Kepler-47 are so dominated by the pri- 
mary star that our analysis does not significantly change 
the approximation of the HZ as compared to a single- 
star model. However, the HZ in Kepler-16 is offset from 
the center of mass of the binary resulting in a planetary 
orbit that alternates in and out of the boundary. The bi- 
nary stars within Kepler-34 and Kepler-35 are similar in 
mass and at such distances from each other that require 
our analysis to determine accurate HZs. The benefit of 
analyzing these Kepler systems, and Kepler systems in 
general, is that they exemplify the diversity we expect 
when exoplanets orbit binary systems. 

To date, these circumbinary planetary systems are 
second only in complexity to the recently discovered 
double-binary system with an interposed exoplanet 
(jSchwamb et al.l I2012T ) . In this instance, an exoplanet 
revolves around an eclipsing binary in a P-type configu- 
ration with a period of ~ 138 days, while a second visual 
binary orbits at ~ 1000 AU. The method that we have de- 
scribed here for determining the HZ boundaries could be 
applied to this quadruple system, as well as any system 
with an n-number of stars. The flux from the n-stars 
would need to be combined and the effective tempera- 
ture determined. However, whether an exoplanet would 
be able to maintain a stable configuration within that 
system would be another matter. We anticipate further 
discoveries by the Kepler mission of new and interesting 
stellar and planetary configurations on which to apply 
our method. 
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